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Abstract 



The quest for life in the Universe is often affected by the free use of extrapolations of 
our phenomenological geocentric knowledge. We point out that the existence of a living 
organism, and a population of organisms, requires the existence of available energy or, more 
precisely, available power per unit volume (§1). This is not a geocentric concept, but a 
principle that belongs to the foundations of thermodynamics. A quest about availability 
in the Universe is justified. We discuss the case in which power comes from mining (§2), 
and from thermal disequilibrium (§3). Thermal disequilibrium may show up in two ways: 
on planets without a star (§4), and on planets where the surface thermal disequilibrium is 
dominated by the incoming photon flux from the nearest star (§6). In the first case we study 
the availability by simulating the structure of the planet with a simple model that contains 
the general features of the problem. For the first case we show that the availability is in 
general very small (§5). In the second case we show that the availability is in general large; 
the order of magnitude depends first of all on the star's temperature and the planet's orbit, 
but is also controlled by the greenhouse gases present on the planet. 



1 Introduction 

We consider the question of whether there is life elsewhere in the Universe, and not 
only on Earth. This question is complicated by the fact that we do not have a defi- 
nition of life as we do for matter, the object of physics: the laws of physics refer to 
something that is the same everywhere and at all times. 

For life we need an ad hoc hypothesis. For instance, we can assume that life is a 
phenomenon unique to the Earth or that life has a uniform chance to exist in the 
Universe, provided certain conditions are satisfied. We take the second hypothesis 
[[l]]. With regards to the conditions, we do have several descriptions of living organ- 
isms on Earth, hence conditions that are nevertheless local and possibly restrictive. 
In fact a certain kind of life, somewhere, could be different in appearance and per- 
formance, adapted to non-terrestrial environments. 

In the absence of a general definition that would resolve these questions, we limit 
ourselves to discuss the most general requirement: available energy. To this purpose, 
we adopt this weak definition of life: a living organism occupies a finite domain, has 
structure, performs according to an unknown purpose, and reproduces itself; in order 
to satisfy these requirements it uses a certain amount of power density (-%)■ That 
a supply of power is necessary, is obvious, upon considering what a living organism 
is not. 

A reversible event does not require power in order to happen; it only requires an 
initial condition, the Hamiltonian does the rest. 

A dissipative structure requires power; examples are the Hadley cells Q in a plane- 
tary atmosphere, which contain ongoing cyclic motions plus chemical reactions, and, 
in the laboratory, the Benard cells ||, the Belousov Zhabotinsky cycle [|j, etc. These 
physical structures have a property: by changing the energy flow across a certain 
boundary, the structure may appear or disappear: there is a bifurcation in a dynam- 
ical dissipative system. 

A living organism is not like this, because it is a control system: its structure is 
related to the operation of controls, that are actuated according to an unknown pur- 
pose. The purpose includes the preservation of the structure itself, plus reproduction, 
plus a certain performance. These concepts exist also in physics and engineering, 
in that branch called control theory. For example a spaceship contains several aux- 
iliary jet engines, that are operated according to a certain strategy of flight. The 
operation of such engines has the effect of changing, at appropriate instants of time 
ti, the initial condition, or phase space configuration at times tj, of the Newtonian 
equations of motion of the spaceship; this behaviour is finalistic rather than deter- 



ministic, even if at all times the spaceship obeys the deterministic laws of gravity. 
Similarly the living organism performs according to a finalistic strategy although at 
all times obeys the laws of physics and chemistry. Both the spaceship and the living 
organism require a power supply. 

In conclusion we make only the hypothesis that the energy requirements for life is 
the same everywhere in the Universe. 

Closely related to this hypothesis is another question that concerns the material 
flow processed by the living organism, once the power is supplied. In fact a living 
organism receives energy and rejects energy, receives chemical elements and rejects 
chemical elements. The chemical elements are extracted and rejected from and into 
the same physical ambient. Such a process can be open, or closed through an am- 
bient cycle connected to the chemical processes occurring in the inside of the living 
organism. This cycle exists on the planet Earth. We do not know much about this 
cycle, therefore it is not clear how its existence can be incorporated as an explicit 
requisite for life everywhere. 

We discuss briefly the absence of the cycle in §2, dedicated to mining, and we limit 
ourselves in the other sections to the evaluation of the available energy. 

Last remark. In the absence of a general definition of life, we have adopted in 
this paper a weak definition that a living organism is finite in space, has structure, 
performs, and reproduces itself. Hence life has to be chemical life, and that is hardly 
conceivable not in connection to a planet. 

An even weaker definition of life is no definition; in this case the search for available 
energy can be pushed to the extreme condition that can be guessed in the asymptotic 
behaviour of the cosmic expansion ||, 0. The available energy, and relative infor- 
mation, come from strategies of mining, all in the vicinity of the cosmic background 
temperature. We call this abstract life. 

It is not clear if an intelligence-preserving transition between chemical life and ab- 
stract life is conceivable in principle. 



2 Energy by mining 

The unit organism occupying a volume V needs a certain power (p. A population of 
n(t) organisms requires a power 

q(t) = n(t)<Pi ([9]=W) 

The population number n(t) is a variable because typical of life is reproduction. The 
time dependence n(t) is determined by the available resources and by the collective 
strategies of evolution of the living population. Mining for energy is a well known 
concept to humankind for about two centuries; that is since the time of mining energy 
from coal and petroleum reservoirs. The population of living organisms functioning 
according to some unknown strategy, extracts from its neighborhood the incoming 
chemical elements a, b, c. . . feeds its own complex state C with both mass and energy 
and ejects the outgoing chemical elements a, /3, 7, . . . to another neighborhood. 

a, b, c I — -(c) - -* g, (3, y 
The energy balance is described by the following figure: 
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Qo and q(t) are expressed in Joule. 

Next we ask which time scales are pertinent to the permanence of the population 
n{t). We wish not to commit ourselves to wild hypothesis on the methodology of 
proliferation of the unknown form of life, the complex state C, and its population 
n{t). The minimum hypothesis is the hypothesis of logistic growth J^J. In this case 
q(t), obeys the logistic equation 

q = aq 9±z±. [g] = J (2.i) 

Qo 

which has the solution 
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The graph of (2.2) is of the kind 
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The collective growth parameter o is the inverse of the time scale r, 

1 

a = — . 

T 

The power used by the population is 

q = n(t) ip; [cp] = W 

and is obtained by derivation with respect to the time of (^ 
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The behaviour of fl2.5| ) is of the kind 




(2.3) 



(2.4) 



(2.5) 



Notice that 
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(2.6) 



Therefore there is a trade-off between mean life r and maximum population n 
From Q and Q: 
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A violation of the rule (|2.7|) can very well be imagined; for some unknown reason the 



population growth may not obey the logistic equation (2.1), which implies inexorably 
the asymptotic exponential disappearance of life. A particular population of living 
organisms may escape the logistic rule and choose survival rather than growth. This 
would be a very interesting life form to discover. 

As a conclusion, we emphasize that the time of survival depends on the strategy of 
the living community itself, is not a property of the ambient in which the hypothetical 
life form develops. 



3 Access to thermal power 

Next we consider an ambient which provides a local flow of power. A local flow of 
power is embedded in a standing, non uniform temperature and pressure distribu- 
tion. For simplicity, we will consider in the foregoing sections only a temperature 
distribution. Two cases can be considered: a non uniform distribution of tempera- 
ture can be a transient, from an initial state to a final state, or can be a steady, non 
equilibrium state. 

A natural realization of the first case is a planet without star during the thermody- 
namical process of cooling off. In fact we may consider the situation whereby after 
several events that brought the planet to its formation and maturity, the structure 
is relatively fixed and the thermal flow develops with time, as in an initial condition 
problem, giving rise to a long transient of cooling. 

A natural realization of the second case is a planet with a star; it receives a flow 
of radiation from the star, and if this flow can be considered constant (for this the 
orbit must be circular) the temperature distribution on the planet's surface, which 
is never uniform, is nevertheless globally stationary in the sense that the total in- 
coming radiation is equal to the total outgoing radiation, and both are constant. 
If the orbit is far from circular, the total incoming flow is no longer constant and 
the temperature distribution can be very complicated. In all cases, any non equilib- 
rium distribution of temperature (and pressure) contains in itself a certain amount 
of "available energy". This concept has been formulated by Gibbs in equilibrium 
thermodynamics ||, and extended by Landau to the formalism of non equilibrium 
thermodynamics H . Here follows a brief review of the concept of Gibbs availability. 
To any situation of thermodynamical disequilibrium there corresponds a well defined 
amount of energy. This is the minimum energy required, in principle, to produce the 
disequilibrium, the uniform equilibrium being the natural state; or, conversely, this 
is the maximum energy which, in principle, could be extracted exploiting the given 
disequilibrium. The words "maximum/minimum in principle" mean that either op- 
eration is done reversibly. Let us see how. 

In equilibrium thermodynamics the concept of Gibbs availability comes by consid- 
ering a closed system in a volume Vq with equilibrium values of pressure and tem- 
perature respectively Po an d Tq; and a portion of volume V < Vq which is out of 
equilibrium with P > Pq, T > Tq. This disequilibrium situation has been created 
reversibly by feeding this amount of energy: 

AW = AQ T ~ T ° +AV(P- P ) (3.1) 



AQ is incoming, as seen from the subsystem V, 

is outgoing, as seen from the system Vo- 
AV is volume reduction, as seen from the subsystem V, 

is volume expansion, as seen from the system Vo- 

Notice that ~ T ° = rj is the Carnot efficiency. 

In the natural return to equilibrium the quantity AW is entirely dissipated, but if 

the disequilibrium situation is undone reversibly, AW could be extracted. In this 

sense the minimum work W becomes the maximum availability. 

In non equilibrium thermodynamics we do not consider two pieces V and Vo, but a 



distribution T(f,t), P(f,t). The generalization of (3A) is intuitive. Rather than a 
quantity of energy AW we consider a power density 

AW 
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The natural extension of (ItT 
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is the gradient operator 
is the fluid velocity 

is the heat flow, 
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It is easy to convince oneself that ( p.2[) is correct. Consider the volume dV = dS"d^ 
During the time dt the increment dQ pertinent to d^ is -H and the Carnot efficiency 




Figure 3.1: 



is r] = — , where dT is the temperature change along d^. Here rj is the usual 
dimensionless efficiency. 



We can write 

W ~ tit dSd£T~ Q V [6 ' 6) 

In(p|) 

VT _, 

ff=—; [ff\ = m 

and (f = d g dt is the common definition of heat flow. In terms of the distribution 
T(f,t), equation fl3.3| ) becomes 

w carnot = fc 0^D_ (3.4) 

We have simply used the definition, valid for heat conduction, 

q = -kVT 

where k is the Fourier constant. 

Remember that if we neglect convective phenomena, the non equilibrium dynamics 

is governed by the Fourier heat equation 

dT 
c v p— = -divq = fcV- VT = kAT 
at 

k is the Fourier constant mentioned above and e = c v pT is the energy equation of 

state. 

For the pressure part, using the same d^, we have, during the time dt 

dPdSd 2 £ = dW; [dW] = J 

d 2 £ is the contraction of d£ (see fig. 3.1), or, equivalently, dSd 2 ^ = d 2 V the contrac- 



tion of dV. There follows 

„ w dW W 

w f lmd = ^ =VP-v; VP • v\ = -^ (3.5) 

dSd£dt ' L J m 3 K ' 



and ( |3.2| ) is justified. 

In the next sections we will limit ourselves on the consideration of yj carnot , equation 

(0). 



4 Planet without star. The transient 

A planet without incoming radiation can belong to a star on a very distant orbit, or 

to a double star on an outer orbit, or can be a planet traveling in some pan- galactic 

orbit. 

For such a planet there is the possibility of containing a non equilibrium global 



transient. In fact, following the idea that planets are formed by accretion |1C], 
we may consider a stage of maturity in which the structure is stabilized and the 
initial temperature is generally a distribution which decreases from the center to the 
periphery. This is a universal feature that has two reasons to occur. 
The pressure must be higher in the deep interior with respect to the surface, due to 
the force of gravity; on the other hand, the two response functions (^f) v and (fy)y 
must be positive, because of the second law. So energy and pressure increase with 
T, therefore the combination of gravity and the second law induces a hot interior. 
The second reason is that the collisions during the process of accretion develop heat, 
and these violent events are followed by a spontaneous creation of layers. In fact 
the various materials present in the planet have in general different heat conduction 
parameters; the good conductors keep transferring heat, the bad conductors remain 
on the surface (T ~ 0), where they cool off but do not transfer much heat and 
therefore constitute the insulating crust. This stratification can be quite complicated, 
although it is a general feature. 

Here we consider the simple model of a core perfect conductor, and an insulating 
crust. The external surface cools off with the Stefan-Boltzmann law 

W 
q = crT ; [q] = — 

var 

q is the heat flow and a ~ 5.67 • 10 -8 m ^. 4 is the Stefan-Boltzmann constant. Ac- 
tually the heat emitted by radiation depends on the nature of the surface, which is 
not necessarily a black solid. If the surface is solid but not black, we may write 

q = e<jT 4 (4.1) 

where e is very close to 1. 

If the surface is gaseous, its molecules may undergo events of inelastic scattering 
with the photons originated by the solid part below. So the gas layer acts as an in- 
sulating greenhouse atmosphere. Nevertheless, in the absence of a star, the surface 
will eventually become cold and the gas phase will not be sustained and we return 
to (4.1) with e ~ 1 valid for solids. 



A general analytic treatment, even assuming spherical symmetry of the temperature 

10 



field, is very complicated. In particular, it requires a guess on the initial condition 
T(r,t = 0); and this has a large uncertainty; secondly it requires a guess on the 
external surface behaviour of VT, and this is again open to uncertainty. 
Our choice is to make the "minimum guess" compatible with the coexistence of two 
problems: what gives the time scale of the cooling process and which is the order of 
magnitude of the Gibbs availability? We anticipate that we will find a relationship 
between availability and lifetime of the transient at the end of §5. 
The model planet we are considering is shown in this figure 




T = T 1 (t) 0<r<R 1 
T = T(r) R l <r<R 2 
T = T 2 (t) r = R 2 

The core contains the thermal energy reservoir at uniform temperature T\(t). The 
insulating layer produces \/T(r) and consequently gives implicitly T 2 (t). The radial 
flow across the shell is 



kr , 2 dT = dl__ 
dr dQ 

dS7 is the solid angle 



const. 



[0]=W 



(4.2) 



We evaluate ^ at r = R 2 
dcf) 



dQ 
so that 

kr' 
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dS(R 2 ) 
becomes 

oTlRl 



aT£R 2 2 



(4.3) 
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The function T(r) is obtained integrating fl4.3|) . We get 

Tir)^- 1 ^^-^ (4.4) 

Notice that T\ and T 2 are at this point unknown. They are related by the following 
algebraic equation, arrived at by putting r = R 2 in Q4.4p 

T 2 = T 1 -aTi (4.5) 

where a is given by 

a = l^lAR [a]= K- 3 (4.6) 

k R\ 

and AR = R2 — R\. The time dependence T\(t) and T 2 (t) is given by the energy 
balance equation: 

-nR\c v pf l (t) = -AirRloT%{t). (4.7) 

This equation simply relates the energy loss from the volume V\ = 3 7r R\ to the 
energy flow across the surface Air R\. Equation ( |4.5| ) gives T\ as a function of T 2 

T x =T 2 + aTt, (4.8) 

putting (^~8|) into (|47|) we get 

J cv p R\ f 2 (1 + 4a Tf) = -i?2 a rl (4.9) 

which is a first order differential equation in T 2 (t) with initial condition 

T 2 (0)=T 20 . 

We may like to take for initial condition the temperature of the core; from ( |4.8| ) we 
may set 



7i(0) 



' 10 



but it is harder to get T 2 from T\. 

Giving our attention to T\ means emphasizing the accretion heating and the total 

mass of the planet. Giving our attention to T 2 means emphasizing the properties of 

the crust, namely its thickness AR and its conductivity. 

Equation (|4.9|) is solved by separation of the variables, T 2 and t. The solution is 



-4a6 1n^- + - (-L - -L I (4.10) 

T 20 3 \T$ Tf ,' 
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where a is given by ([O]) and b is given by 



SaRj' 



[b] = K 3 • s 



(4.11) 



Equation ( |4 . 1 Q|) implicitly gives T 2 (t;T 2 o), or, in brief, T 2 (t); and using ($L 



we 



have Ti(t;r 2 o), or Ti(t). With this knowledge, equation (4^4) is now completely 
determined and we can write 



T(r.l) - '/',(/:'/:,,,) - T <rTi(t;T 20 )R 2 2 "—^ 



1 



(4.12) 



Consider the mean life of T 2 (t) coming from (4.10), namely the time r such that 

(4.13) 



T 2 (t) = e 1 T 20 



We get 



T = 4ab + \(e 3 -l)b 7 ± r 



n + t 2 



(4.14) 



The first addend n contains only the planet's parameters, while r 2 contains also the 
initial condition. This suggests that the differential equation (|4.9| ) has two regimes. 
In fact, consider the domains of T 2 such that 



T 2 > 



4a 



- ^" 



a) 



6) 



In the range ( 4.15| a), equation ( pL9| ) is approximated by 



-ac v pR(T 2 = -R 2 a T 2 



The solution of ( 4.16| ) is 

T 2 (t)=T 2 (0)e~n 



(4.15) 



(4.16) 



(4.17) 



with 



Ah 4 °"P A » ^1 

Tl = 4a6= 3— Ai? i^ 



The value 
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fclfr 



4aR 2 AR 



(4.18) 



(4.19) 
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is relevant also for equation ( [4 . 8| ) : in fact in the range ( 4.15 b) the relationship 



between T\ and T 2 is linear, while in the range ( 4.15| a) the relationship is dominated 
by the behaviour 



Ti~aT 2 4 ; 




1 the range (4.15|b) equation 


1 • 
- c v p R{ T 2 = 


-RJaTi 


ith solution 




\ J 20 


3 \ _ * 



is approximated by 



(4.20) 



(4.21) 



(4.22) 



-'-s^-'ifsa < 423 > 



and 

1 (J> i\ c vP R l 1 

? 2 T 3 

l 2 J 20 

Remember that in this model the planet is never at uniform temperature for R\ < 

r < i?2, but the assumed initial distribution is T\ S> T2, so in general it is the 

exponential solution fl4.17| ) that interests us. 

We show in fig. [4.1| the function T\ = T\{T2) or T 2 = T 2 (Ti) for a given choice 

of the parameters k and AR; a is a fundamental constant. The region T 2 < T 2 is 

magnified. 

Notice that T\ and T 2 have different scales and the linear region T\ = T 2 appears 
flat. For T 2 < T 2 we also have that V r T is small. 

We now proceed in this section adopting an approximation which is not physically 
restrictive and simplifies the calculations. We assume that the crust AR of the planet 
is small with respect to both R\ and R2 

R2 t , 

-L „ „ 1 (4.24) 

Using (|4.24|) , equation ([4.12| ) becomes linear in r: 

T(r,t) = Tx(t) -l<rT$(t) (r - R x ) (4.25) 
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Figure 4.1: The relationship T\ = Ti(T 2 ) and T 2 = T 2 (Ti). In this figure we have taken k 

0.3 W m _1 IT 1 , AR = 20 km. It follows from fa.lfy T 2 * ~ AK. 



T x = T 2 + a'T£, 



(4.27) 
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(4.28) 
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(4.29) 



From now on we adopt the linear approximation, namely we use the above expres- 
sions. However we do not continue to carry the primes. 



15 



5 Thermal availability in the shell 

We have now the ingredients necessary to evaluate the Gibbs availability in the shell. 
With reference to §3, we use the Carnot availability 

(VT) 2 
w = k^- (5.1) 

and we have for our model planet 

T(r,t)=T 1 (t)-^a(r-R l )T 2 4 (t) (5.2) 



(5.3) 

(5.4) 
We see that w(r, t) increases with r in the interval R\ < r < R 2 and we obtain from 

for any given t (5-5) 



VT = — = 

dr 


-\°T$(t) 


so that 




w(r,t) = — 





(5.4) this ratio 

w(Ri) _T 2 

w(R 2 ) T x ' 

Let us discuss now the trade-off between availability and thermal mean-life of the 
host planet. In the foregoing we consider the regime T 2 > T 2 * ; the mean life r ~ t 2 
no longer appears, so we will use the symbol r for t\ in what follows. 
We consider 

w(r,t) = T{r ^ (5-6) 

in the regime T 2 > T 2 *, that means that 

Ii~(|)\ (5.7) 

so that (15.61) becomes 



i^W 



"''-■') = r(r,«) (5 ' 8> 

Using the linear approximation 
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we now have 

w(r,t) = k^ fy (5.9) 

At a given T2 and for fixed r between Ri and R2, the dependence of w on AR 
satisfies the following scaling law. 

w(a AR) = -^ w(AR) (5.10) 

or 

Always in the regime T<i > T£, we have 

r = ^R 1 AR (5.11) 

At a fixed r we have 

T(aAR) = ar(AR). (5.12) 

The comparison between ( |5.10| ) and ( |5.12| ) gives the following result: for a given 
supply of thermal energy stored in the planet, | ix c v p Rf T\ , a larger availability in 
the shell implies a thinner shell and therefore a shorter mean-life of the stored energy. 
As far as we know from the properties of terrestrial life, living organisms require a 
rather narrow range of ambient temperature. For example they need liquid water. 
The phase diagram of water is known and gives the domain T, P where we have the 
liquid phase. To quantify these ideas, we define a best temperature Tf, and we set 

T fe = 300K 

Fixing Tfe poses a problem because both T and w are function of r and t. We wish 
to show that the imposition 

T(r,t)=T b (5.13) 

has implications which can be explicitly calculated. First of all ( |5.13| ) is an algebraic 
equation in r which defines a function 

r b = r b (t) (5.14) 

and, putting ( 5.13| ) into ( ^.6[ ), defines also 



w b = w b (t) (5.15) 

In particular we will show that 

r b (t) < r b (0) for t > (5.16) 
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W b {t) < 1^(0) 



for t > 



(5.17) 



The inequality ( |5.16 ) implies that the conditions for life must migrate from a loca- 
tion near the surface toward the interior, and correspondingly ( |5.17| ) says that the 
availability decreases with the migration. Equation ( |5.14 ) can be written explicitly 
in the linear approximation 

(Ti(t)-T b ) 



n(t) = Ri + 



m 



(5.18) 



We see that (obviously we consider T 10 > T b ) T\ reaches the value T b when t reaches 
a certain value t, namely 



n -> Rr, 



r, 



for 



t -+t 



In other words for t — > i the location of life reaches the bottom of the crust. From 
that moment on, T\(t) keeps decreasing and the requirement T — T b can no longer 
be maintained. Also the availability can be written explicitly 



w b 



a -T$(t) 



(5.19) 



and we see at once that w b decreases with t. 

The question is whether t is large or not. We have defined the time t with the 

equation 



T 1 (t)=T b 
Using 

T 1 (t)~aTi(t) = aTi e- 
we obtain 



4± 



t ~ - In — -^ 
4 \ T b 



4 V^fe 



(5.20) 



(5.21) 



(5.22) 



Notice that 



t = 



if 



•10 



T b 



(5.23) 



l > r 
< 



if 



iHlflM 



In l'^ I < ! 



(5.24) 
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We can write the time dependence of r b and w b explicitly, showing how it is affected 

by the model parameters. 

We are in the T 2 > T% regime, so that 

r 2 (t)~r 20 e-7 (5.25) 

with 

t~ -^AR^ (5.26) 



Using (fT25D in fl5lg| ) and (|5l9D we get (for < t < t): 



and 



r&( 


= 


-Ri 


ka 
+ — 
a 






-Rx 


ka 
+ — 
a 


«',, 


(*) = 


a 2 
k 


J 20 

T b 6 



k T b 4 t 

a J 20 



(5.27) 



(5.28) 



It turns out, by inspection of ( 5.27] ) and ( |5.2Sj ) that the crucial parameter is r. 



Knowing the value of r, we can estimate the migration length r b (0) — r b (t) = 5r(t) 
and the availability decrease w b (0) — Wb(t) = 5w(t) after a certain time t. From 
(|5 . 2 7| ) we have 

5r b (t) = * £- (e^ - l) (5.29) 



while from (5.28) we have 



Note that the linearization of T{r, t) and the limitation to T 2 > T 2 * are reasonable 
simplifications. Obviously r b (t) and w b {t) could be calculated numerically choosing 
the parameters of the model and the initial condition T 2 o at will, and without the 
above approximations. 

We conclude this section by giving a numerical estimate of r b , w b , 5r b , 5w b and 
t for a certain set of the planetary parameters. We take a planet that vaguely re- 
sembles the Earth. 

Material parameters 
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c v 5 • 10 3 J kg^ 1 KT 1 Iron 
p 10 4 kg m~ 3 Iron 

k 0.3 W m" 1 K" 1 Sand soil 

Geometry 

R x 6.28 -10 6 m 

R 2 6.3 -10 6 m 

AR = R 2 - Ri 20 km 

a is the Stefan-Boltzmann constant ~ 5.67 • 10~ 8 W mT 2 K -4 , Tj, = 300 K by 

our choice. 

Initial condition [0]: 

W 

q ~ 0.05^ 

Note: The Earth is about 10 9 years old. Here we take the present age as the initial 
time. Alternatively we could take a higher qo and correspondingly higher T\q and 
begin our clock much earlier. The reader can play with these numbers. 
Since 

qo = a T 20 ; we have T 2 q = ( — J ~ 30.6 K; 

and 

aTL = 3314K 



a = ^Ai2 = 3.78-1(T 3 K- 3 

k 

We calculate rb — R\ for t = 

r 6 (0) -R 1 = - a (l ) ~ 18.2 • 10 3 m (5.31) 

0- V T wJ 

which means r^ is ~ 2 km below the surface at r = R 2 , and Wb for t = is 

— T 8 w 

Wb = J*-* „ 2.75 • 10~ 5 -^. (5.32) 

T b m d 

Concerning the time dependence we calculate r and get 

T ^ - ( bL£-R 1 /\R~ 8.85 • 10 11 yr. (5.33) 

3 k 



Given r, using equation ( |5.22 ) we calculate i 



T • f T w \ . . c q mil 



i = - hi I -^- ) ~5.3-10 1J -yr~0.6r. (5.34) 
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What are the quantities Sr^ and SiVb after t = 10 6 yr? We get from ( 5.29| ) and ( |5.30| ) 

l) ~ 0.008 m (5.35) 

* 20 v ' 

From (|5.28|) we have 
<^(10 6 yr) 



^(10 6 yr) = ^^(e 4 7 



a 2 n 



20 

k T b 



1 



2.5- 10 -10 



W 



m" 



(5.36) 



The faintness of the continuous availability evaluated in the above numerical exam- 
ple is an indication of a general property. In fact we cannot make wild hypotheses 
with regards to T20, Ai? and the core equation of state. In particular, AR governs 
the trade off between the planet's thermal life, r, and the availability w that it can 
house. It is therefore obvious that w is in general small. 

We may consider the possibility of surface availability, rather than volume availabil- 
ity. This implies a discontinuity of T, and VT, and this may exist in the presence 
of a small scale structure on the shell of the planet, for instance some form of hot 
streams coming from the core and protunding into the cooler crust. In this case 
we may have an interface between the stream at T\ and the ambient at T%. Such a 
discontinuity produces an availability 



w = g-lY^— ^W-£o) 

where we make reference to the following figure: 



(5.37) 




The quantity <f is the heat flux, £ is the normal vector, the straight line £ being per- 
pendicular to the surface element dS, T\ is the stream temperature, T2 the ambient 
temperature. 



The justification for considering the surface availability (5.37) is that temperature 
discontinuities exist on Earth. We have various examples of hot springs, both at the 
interface solid-atmosphere and solid-ocean. We may reasonably think that this kind 
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of structure may exist also on other planets. 

On Earth, with such discontinuities there sometimes corresponds a local commu- 
nity of living organisms, an opportunistic niche. Can we make the hypothesis that 
such niches exist on planets without a star? This is not an innocent hypothesis, 
because it implies the preliminary verification of a very strong statement, namely 
that the existence of the terrestrial niches attached to discontinuities is disjoint from 
the existence of the global living community (which for our planet is a photon driven 
system) both in an ecological and evolutionary sense. In the absence of a proof of 
this statement, the consideration of life niches belonging to discontinuities on planets 
without star is not a tenable hypothesis, but a belief. 
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6 Planet with star. The global steady state 

A planet that receives radiation from a star has on its surface a thermal distribution 
much different from that of a planet without star. What we have discussed in §4 
is a global transient with spherical symmetry. The interesting dynamics is in the 
insulating shell. 

If a planet is sufficiently close to its star, the radiation flow impinging on the planet 
can be stronger than the energy flow originated by the planet's core, and the planet's 
surface dynamics becomes dominated by the external input. 

Note. Numbers for the Earth: power flow from the core go = 0.05 -%, solar constant 
^ = 1350-^. 

The external input is no longer symmetric in the reference system of the planet. 
On the other hand, if the total flow investing the planet is constant, the surface 
of the planet can be in a global steady state. The first question is therefore with 
regards to the astronomical conditions necessary to have a global steady state. 
In the solar system the planetary orbits are in good approximation elliptic and also 
with small eccentricity. A good astronomical question is whether this architecture 
is unique to the Sun and its planets or is a general Newtonian property. We leave 
the deep and difficult problem of the stability and structure of planetary systems to 
the astronomers. We remark, however, that in the absence of near circular orbits 
the radiation flux hitting a planet may be a chaotic (non periodic) function of the 
time [O ], and the thermal distribution, created on the surface of the planet, will be 
something which is neither a transient nor a steady state. 

Consider now a near circular orbit. The next question concerns the temperature 
of the star T s , with radius R s and the distance star-planet d. In fact from these 
three numbers follows the average surface temperature of a planet with a black body 
surface and without an atmosphere. The average temperature on the circular orbit 
of radius d is calculated from the energy balance. In this formula (j) m is the portion of 
the isotropic total radiation emitted by the star and intercepted by the cross section 
of the planet ir R 2 . 4> out is the total energy radiated by the planet. 

r>2 

cj) in = 7TR 2 aT 4fl£ = ^out = 47T R 2 a T* v (6.1) 

From the balance (|6.1|) we get T av : 



1 / 7? 

Now T av <C T s because the radiation flow originated on the star is diluted at the 



distance d; but the spectroscopic temperature remains T s 12]. 
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In turn, if we look for long stellar lifetimes inside the main sequence, T s , R s and t s 
are related to the stellar mass. If we require a long duration of stability, for instance 
t s ~ 10 9 yr, the values of T s and R s are not wildly different from those of the Sun, 
so the consideration of T s in the order of few thousand K, makes sense. 
A global steady state takes place if the global balance ([O]) is satisfied with a surface 
temperature T(6, (p, t) which depends on position and time. The detailed properties 
of the global steady state are related to the total amount of greenhouse gases. In 
fact 4> out can differ from (f) m during a transient. This transient can only be a period 
of increase or decrease of greenhouse molecules in the planetary atmosphere. Now, 
except for external causes, primarily planetary collisions, there is no good reason 
to believe that the total amount of greenhouse gas must oscillate. It follows that 
a global non equilibrium steady state is the rule. If the amount of greenhouse gas 
is constant, the effective average surface temperature Tal will be different from 
the average temperature T av , calculated above, that takes into account only the 
astronomical parameters. We can in general say that 

rpeff _ rp 

A av — 1 ± av 

with 7 > 1. The value of 7 depends on the history of the planet's formation, and 
the actual amount of greenhouse gases existing at a certain stage of its evolution. 
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7 Photon availability 

At present we do not have any ideas about the possibility of gaseous life. Life is 
structured, shaped; gases are, by definition, interactionless and shapeless. So we 
disregard it here and consider the photons that have passed through the atmosphere 
and hit a liquid or solid layer of molecules that can thermalize. We may assume 
that only a few collisions are sufficient to thermalize the radiation flow, or, in other 
words, the thermalization is a surface effect. The surface temperature changes with 
time and the latitude on the planet: day and night, and the seasons. This is a 
very important fact for the understanding of the complexity of photon-powered life. 
T(f,t) oscillates around the value Tal , the average value that takes into account 
the greenhouse effect. 

The photon availability vf 1 is rather different from the conduction availability, yj carnat } 
studied in the preceding sections. In fact, considering that the layer involved in the 
thermalization of the photons is very thin, we may write 



uP{r) = (p f] 5{r — R) 



(7.1) 



{8 is the Dirac function, the symbol 7 stands for photon). 



w 



w 

III 3 ' 



[5(r)} 



m 



-1 



In (7.1) 



f? 2 



m<p,t) 



(7.2) 



where f(6, ip, t) is a kinematic factor depending on latitude and longitude. The factor 
f(9,(p,t) takes into account the relative motion between star and planet. See ref. 
|12|| for details. The efficiency -y- is given by 



V 



AT 
~T~ 



T(R) 



(7.3) 



The efficiency r\ is an astronomical property rather than a planetary one. In fact, 
using for T(R) the value T av (|J) 



T(R) 



we get 



T — — i Rs T 
av ~ v/2 V d s 



1 

V2 




We see immediately that 




R) 



(7.4) 



(7.5) 



(7.6) 
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Equation ( |7.6|) says that the photon availability is an astronomical property, and the 
planet radius appears uniquely in the delta function 5(r — R). The photon availability 
has several very important properties. 

- T s is high, and is in a range that may involve important chemical reactions 

- The power input w" 1 works mainly on the surface, 

- There is in general a high value of rj, namely r/ ~ 1, or, in other words T s and T av 

are quite different from each other. 

As citizens of the planet Earth, we are the best observers of life powered by photon 
availability. We see that 

- An organism powered by photons is richly structured. This is a direct consequence 

of the fact that the living organism must develop itself behind the receiving 
surface: it must be organized in a structure. 

- The photon-driven chemical reactions constituting the input to the complex state 

C, can be, in principle, the alternative to mining of chemical components. This 
alternative is indeed realized by Nature. In fact it happens that the complex 
living state C, which does not need mining in the input, does not need a 
reservoir for the outgoing chemical states either. 

We have described the metabolism based on mining with the graph (§2) 
a, b, c | — -(c) - -* g, (3, y 

The complex state C of the terrestrial ecosystem belongs to a different graph: 




a, b, c oc, p\ y 

'""■■ — ^ -■■* 

On the surface of the planet Earth the outgoing components a, /3, 7 circulate in the 
physical ambient and perform a cycle reproducing exactly a, b, c in the amount that 
is required by C. It is uniquely by closing the cycle a, b, c — > a, /3, 7 that life on a 
planet can be eternal, where eternal means over the lifespan of the host star. 



26 



Notice that the cycle has an upper limit to its total power and this is the integral of 
the photon availability on the planet's surface 

$ mm ~ / q(e, V, t) Ts - T ( 0, ^ ^ dS; [$ m ^]=W (7.7) 

Js T s 

T(8, <p,t) is the temperature field on the planet's surface. Obviously the total avail- 
ability cj> max is utilized by the population of the living organism only in part. The 
other part goes into the physical channel, mainly the non-organic water cycle and 
thermal circulation. We know that on Earth about 10 -2 of the total incoming flux 
goes into the life channel [O], This ratio comes from observation, not from theory. 
But the important fact is that life based on stellar radiation and organized accord- 
ing to the above cycle has a physical upper limit to its growth, and consequently its 
duration is not population dependent, as it happens with life from mining (§2) but 
star dependent. 
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8 Conclusion 

Comparing the numerical discussion of the thermal availability in §5 and the photon 

availability in §7, it appears that the physical chances for life on a planet without a 

star are rather miserable in comparison to a planet in a near circular orbit around a 

star. 

In both cases, the time scales for the endurance of a hot planetary core, or the burning 

of the star, are sufficiently long; that is not the limiting effect for the permanence 

of life. It is the material cycle (§1) that matters. It is conceivable that when the 

available energy is very low, the strategy of life goes into the direction of mining, for 

both chemical elements and energy. If this is the choice, life becomes a transient. 

It may exist by bursts, and chances of communication with other bursts in other 

places in the Universe may become negligible. Life, as developed on Earth functions 

according to a global cycle, and therefore the duration of life on the planet can be 

equal to the duration of its host star (§7). 

The understanding of the life-ambient cycle is a much harder problem than the mere 

evaluation of the availability. It is reasonable to think that the existence of this cycle, 

and the corresponding permanence of life, may pose other very severe conditions for 

the possibility of life in the Universe. 

At present we are not able to formulate these conditions: we are only able to witness 

that they exist on Earth. 
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